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SYSTEM AND METHOD FOR HIERARCHICAL ANALYSIS OF CONTRAST 
ENHANCED MEDICAL IMAGING INFORMATION 

TECHNICAL FIELD 

The present disclosure relates generally to medical imaging techniques. 
More particularly, the present disclosure describes systems and methods for 
hierarchically analyzing and/or evaluating dynamic imaging properties 
corresponding to contrast enhanced medical imaging information to aid differential 
lesion diagnosis procedures. 

BACKGROUND 

Medical imaging technologies can provide detailed information useful for 
differentiating, diagnosing, or monitoring the condition, structure, and/or extent of 
various types of tissue within a patient's body. In general, medical imaging 
technologies detect and record manners in which tissues respond in the presence 
of applied signals and/or injected or ingested substances, and generate visual 
representations indicative of such responses. 

A variety of medical imaging technologies exist, including Computed 
Tomography (CT), Positron Emission Tomography (PET), Single Photon Emission 
Computed Tomography (SPECT), and Magnetic Resonance Imaging (MRI). 
Particular imaging techniques, for example, certain MRI techniques, may scan a 
volume of tissue within a region of anatomical interest. Scan information or data 
corresponding to an anatomical volume under consideration may be transformed 
into or reconstructed as a series of planar images or image "slices." For example, 
data generated during a breast MRI scan may be reconstructed as a set of 40 or 
more individual image slices. Any given image slice comprises an array of volume 
elements or voxels, where each voxel corresponds to an imaging signal intensity 
within an incremental volume that may be defined in accordance with x, y, and z 
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axes. The z axis commonly corresponds to a distance increment between image 
slices, that is, image slice thickness. 

Any given medical imaging technology may be particularly well suited for 
differentiating between specific types of tissues. A contrast agent administered to 
the patient may selectively enhance or affect the imaging properties of particular 
tissue types to facilitate improved tissue differentiation. For example, MRI may 
excel at distinguishing between various types of soft tissue, such as malignant 
and/or benign breast tumors or lesions that are contrast enhanced relative to 
healthy breast tissue in the presence of Gadolinium DPTA or another contrast 
agent. 

Medical imaging techniques may generate or obtain imaging data 
corresponding to a given anatomical region at different times or sequentially 
through time to facilitate detection of changes within the anatomical region from 
one scan to another. Temporally varying or dynamic tissue dependent contrast 
agent uptake properties may facilitate accurate identification of particular tissue 
types. For example, in breast tissue, healthy or normal tissue generally exhibits 
different contrast agent uptake behavior over time than tumorous tissue. 
Moreover, malignant lesions generally exhibit different contrast agent uptake 
behavior than benign lesions ("Measurement and visualization of physiological 
parameters in contrast-enhanced breast magnetic resonance imaging," Paul A. 
Armitage et al., Medical Imaging Understanding and Analysis, July 2001, 
University of Birmingham). 

In general, at any particular time, the intensity of an imaging signal 
associated with any particular voxel depends upon the types of tissues within an 
anatomical region corresponding to the voxel; the presence or absence of a 
contrast agent in such tissues; and the temporal manners in which such tissues 
respond following contrast agent administration. In several types of breast MRI 
situations, normal or healthy tissue exhibits a background signal intensity in the 
absence of a contrast agent, while abnormal or tumorous tissue exhibits a low or 
reduced signal intensity relative to the background intensity. Thus, prior to 
contrast agent administration, abnormal tissue typically appears darker than 
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normal tissue. In the presence of a contrast agent, lesions or certain types of 
abnormal tissue typically exhibit a time-dependent enhancement of imaging signal 
intensity relative to the background intensity. 

In general, a lesion will exhibit one of three types of time dependent 
contrast agent uptake behavior. Within imaging data corresponding to a time 
series of scans, each type of contrast agent uptake behavior manifests as a 
corresponding type of dynamic imaging signal intensity profile or curve. Each type 
of dynamic intensity curve probabilistically corresponds to whether the lesion is 
benign or malignant. 

Figure 1 is a graph 100 that generally illustrates a first, a second, and a 
third type of dynamic imaging signal intensity or relative enhancement curve 110, 
120, 130 that may be obtained for a given region of interest (ROI) (i.e., a set of 
voxels corresponding to an anatomical region) that encompasses one or more 
portions of a lesion. In the graph 100, the horizontal axis corresponds to time, 
while the vertical axis corresponds to imaging signal intensity or the extent to 
which imaging signal intensity is enhanced relative to an initial or precontrast 

imaging signal intensity. 

Prior to contrast agent administration, a precontrast scan is initiated or 
performed at a first scan time to. For any given voxel, the precontrast scan 
establishes a precontrast imaging signal intensity and/or a reference relative 
enhancement value that may be represented as a variable S 0 . Contrast agent 
administration occurs some time after to at a time tc, for example, one minute after 
t 0 . Essentially immediately after contrast agent administration, the level of 
contrast agent within vasculature associated with a lesion increases. Imaging 
signal intensity or relative imaging signal enhancement associated with the lesion 
correspondingly increases, typically at a rapid rate during an initial time interval. 
Depending upon lesion characteristics, tissue dependent contrast agent kinetics 
may subsequently give rise to in an imaging signal intensity curve that (1) 
continues to increase or enhance; (2) reaches a peak level of enhancement and 
then levels off or plateaus in an abrupt or generally abrupt manner; or (3) reaches 
a peak level of enhancement and subsequently declines to lower or reduced 



-3- 



levels of enhancement. Each of the aforementioned curve types is respectively 
referred to as (1) a continued, steady, or persistent enhancement curve 110; (2) a 
plateau curve 120; and (3) a washout curve 130. 

A first postcontrast scan is performed at a first postcontrast scan time ti, 
which generally corresponds to a time at or near which a peak in imaging signal 
intensity and/or relative enhancement Si would be expected in the context of a 
plateau or washout curve 120, 130. The first postcontrast scan time ti may be, for 
example, one minute after contrast agent administration. For ease of 
understanding, Si is shown in Figure 1 as having an identical value for each curve 
type 110, 120, 130. More generally, each type of curve 110, 120, 130 may 
correspond to a unique or distinct Si value. 

A second postcontrast scan is performed at a second postcontrast scan 
time t 2 , thereby capturing or acquiring for each voxel another corresponding 
imaging signal intensity and/or relative enhancement value S 2 . The interval 
between ti and t 2 is sufficiently large to improve or maximize a likelihood that 
imaging signal intensity differences corresponding to times ti and t 2 can facilitate 
categorization of an imaging signal intensity curve as a persistent enhancement, 
plateau, or washout curve 110, 120, 130. The time between ti and t 2 may be, for 
example, approximately 4 minutes. 

Dynamic imaging signal intensity or relative enhancement curves are 
typically numerically characterized in accordance with two parameters, namely, a 
percent enhancement (PE) value and a signal enhancement ratio (SER). For a 
given voxel, the PE value is defined as the difference between the first 
postcontrast imaging signal intensity Si and the precontrast signal intensity S 0 , 
normalized relative to the precontrast signal intensity S 0 . The SER for any 
particular voxel may be defined as the difference between the first postcontrast 
imaging signal intensity Si and the precontrast signal intensity S 0 , normalized 
relative to the difference between the second postcontrast imaging signal intensity 
S 2 and the precontrast signal intensity S 0 . For persistent enhancement, plateau, 
and washout curves 110, 120, 130, the SER will have a value that is less than 1 .0, 
equal to 1 .0, and greater than 1 .0, respectively. 
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[0013] PE values exhibit a correspondence to lesion type. In general, a higher PE 

value may suggest a higher probability that a lesion is malignant. Some existing 
systems and/or methods for analyzing dynamic MRI data identify a set of voxels 
corresponding to a highest PE value or a highest intensity as a malignant lesion. 
Other existing systems and/or methods may identify (1) a set of voxels 
corresponding to a PE value above a first reference value as a malignant lesion; 
(2) a set of voxels corresponding to a PE value below the first reference value and 
above a second reference value as an indeterminate type of lesion; and (3) a set 
of voxels corresponding to a PE value below the second reference value as a 
benign lesion ("Dynamic Breast MR Imaging: Are Signal Intensity Time Course 
Data Useful for Differential Diagnosis of Enhancing Lesions?", Christiane Kuhl et 
al., Radiology, April 1999). 

[0014] Curve type 110, 120, 130 also exhibits a correspondence to lesion type. In 

particular, a washout curve 130 is strongly indicative of a malignant lesion or 
weakly indicative of a nonmalignant or benign lesion. A plateau curve 120 may be 
somewhat more indicative of a malignant lesion than a benign lesion. Finally, a 
persistent enhancement 110 curve is strongly indicative of a benign lesion or 
weakly indicative of a malignant lesion. Certain methods for analyzing dynamic 
MRI data identify a curve shape of washout 130 or plateau 120 as malignant, and 
a persistent enhancement 110 curve as benign (ibid). 

[0015] Unfortunately, systems and/or methods that analyze contrast enhanced 

medical imaging data in manners described above fail to adequately increase or 
maximize diagnostic sensitivity, specificity, and/or accuracy. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0016] Figure 1 is a graph that generally illustrates a persistent enhancement, a 

plateau, and a washout type of dynamic imaging signal intensity profile or curve 
that may be obtained for a given region of interest that encompasses one or more 
portions of a lesion. 
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[0017] Figure 2 is a block diagram of a system for hierarchical analysis of contrast 

enhanced medical imaging information according to an embodiment of the 
invention. 

[0018] Figure 3 is a flowchart of a procedure for hierarchical analysis of contrast 

enhanced medical imaging information according to an embodiment of the 
invention. 

[0019] Figure 4 is a flowchart of a procedure for identifying or generating one or 

more volumes of interest (VOIs) according to an embodiment of the invention. 

DETAILED DESCRIPTION 

[0020] The present disclosure describes systems and/or methods for hierarchical 

analysis, evaluation, characterization, and/or processing of contrast enhanced 
medical imaging information and/or data. Such systems and/or methods may 
facilitate Computer Aided Detection and/or Computer Aided Diagnosis (CAD) of 
particular types of medical conditions, tissues, and/or tissue abnormalities such as 
malignancies. Portions of the following description detail manners in which certain 
embodiments of the present invention may be applied in an MRI context, 
particularly MRI imaging of breast tissue. Notwithstanding, various embodiments 
of systems and/or methods in accordance with the present invention may be 
applicable to essentially any type of dynamic medical imaging technology and/or 
technique that utilizes a contrast agent. 

[0021] One or more portions of various systems and/or methods described herein 

may be applicable to two and/or three spatial dimensions. For ease of 
understanding, in the context of the description that follows, a volume of interest 
(VOI) may be defined or construed to include, encompass, and/or equivocate to a 
region of interest (ROI). Similarly, a voxel may be defined or construed to include, 
encompass, and/or equivocate to a pixel. 

[0022] Systems and/or methods in accordance with particular embodiments of the 

invention may analyze a VOI in a selective hierarchical sequence that may include 
a washout analysis procedure, a plateau analysis procedure, and/or a persistent 
enhancement analysis procedure. A VOI may comprise at least one, and typically 
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multiple, candidate voxel sets. The washout, plateau, and/or persistent 
enhancement analysis procedures may respectively determine whether imaging 
information associated with one or more candidate voxel sets under consideration 
exhibits dynamic properties corresponding to or indicative of washout, plateau, 
and/or persistent enhancement behavior. 

A washout, plateau, and/or persistent enhancement procedure may use 
particular values and/or value ranges of one or more types of parameters 
indicative of dynamic imaging signal behavior as curve type indication thresholds. 
Such procedures may evaluate acquired and/or generated imaging signal data 
relative to such thresholds to facilitate identification, characterization, recognition, 
and/or determination of washout, plateau, and/or persistent enhancement 
behavior. 

In order to improve, increase, or maximize diagnostic sensitivity, specificity, 
and/or accuracy, a washout analysis procedure may be initiated, performed, 
and/or executed first, selectively followed by a plateau analysis procedure, 
selectively followed by a persistent enhancement analysis procedure. If any given 
washout, plateau, and/or persistent enhancement analysis procedure determines 
that enhancement behavior associated with one or more candidate voxel sets 
under consideration matches or corresponds to a type of dynamic behavior to 
which the analysis procedure is directed, the analysis procedure may identify as a 
malignancy or a likely malignancy a candidate voxel set that exhibits, essentially 
exhibits, approximately exhibits, or is expected to exhibit a strongest, most 
significant, best, most likely, or most adequate correspondence to malignant 
behavior. In certain embodiments, a voxel set corresponding to a highest percent 
enhancement (PE) may be identified as malignant or a likely malignancy. 
Alternatively or additionally, a voxel set corresponding to a steepest, strongest, 
and/or fastest imaging signal intensity decline may be identified as malignant or a 
likely malignancy. If a washout or plateau analysis procedure has respectively 
identified a likely malignancy corresponding a given VOI, a subsequent plateau or 
persistent enhancement analysis procedure directed to this VOI may be avoided 
or skipped. 
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Figure 2 is a block diagram of a system 200 for hierarchical analysis of 
contrast enhanced medical imaging information according to an embodiment of 
the invention. The system 200 may comprise a medical imaging system 210, a 
data storage unit 220, and a computer system 240 configured to perform one or 
more types of hierarchical imaging information analysis procedures according to 
particular embodiments of the invention. In one embodiment, each element 210, 
220, 240 of the system 200 is coupled to a computer network 290. The medical 
imaging system 210 may comprise an MRI or other type of imaging system. The 
data storage unit 220 may comprise one or more types of hard disk, optical, 
and/or magneto-optical drives, and may possibly comprise a Network Attached 
Storage (NAS) device. The data storage unit 220 may receive, store, and/or 
transfer imaging data as well as other information. 

The computer system 240 may comprise one or more portions of a medical 
image analysis platform. The computer system 240 may include a processing unit 
and a memory, and may further include one or more of a disk drive and/or other 
data storage devices (e.g., optical and/or magneto-optical data storage devices, 
tape drives, flash memory based drives, etc.), an input device, and an output 
device. For the sake of clarity, these conventional and well known computer 
components are not illustrated in Figure 2. The memory, the disk drive, and/or 
other data storage devices may comprise one or more portions of computer 
readable media that store program instructions and possibly data for performing 
one or more hierarchical imaging information analysis procedures and/or 
operations associated therewith in accordance with particular embodiments of the 
inveniion. Depending upon embodiment details, the network 290 may comprise 
one or more local or private networks such as a Local Area Network (LAN) and/or 
one or more public networks such as the Internet. 

Figure 3 is a flowchart of a procedure 300 for hierarchical analysis of 
contrast enhanced medical imaging information according to an embodiment of 
the invention. In one embodiment, a hierarchical analysis procedure 300 includes 
an acquisition procedure 302 that involves acquiring, generating, retrieving, 
receiving, and/or obtaining imaging information, signals, and/or data 
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corresponding to a set or series of medical images. Depending upon embodiment 
details, imaging information may correspond to breast images (e.g., breast MRI 
scans) or other types of images. The acquisition procedure 302 may involve or be 
directed toward at least one precontrast and at least two postcontrast image 
volumes to facilitate numerical characterization of dynamic or time dependent 
imaging signal behavior. As further described below, such numerical 
characterization may involve determination and/or analysis of signal enhancement 
ratios (SERs) based upon time dependent voxel intensities S 0 , Si, and S 2 . In 
certain embodiments, an acquisition procedure 302 may consider more than two 
postcontrast image slices to facilitate analysis of time dependent voxel intensities 
across additional levels of temporal granularity. 
[0028] The hierarchical analysis procedure 300 may further include a VOI 

identification procedure 304 that involves identifying or determining one or more 
VOIs for consideration within the imaging data, where each VOI comprises at least 
one candidate voxel set. In certain embodiments, elements within a candidate 
voxel set correspond to voxels that have a potentially increased, increased, and/or 
significant likelihood of corresponding to abnormal or malignant tissue. 
Depending upon embodiment details, a candidate voxel set may include one or 
more individual voxels and/or one or more combined voxels, where a combined 
voxel may comprise a plurality of voxels that are averaged, mixed, and/or blended. 
A combined voxel may be based upon or generated through one or more 
operations upon imaging signal information within or corresponding to a 3 x 3 x 1 
voxel group; a 9 x 9 x 1 voxel group; a 3 x 3 x 3 voxel group; and/or an otherwise 
dimensioned voxel group. One or more candidate voxel sets may be identified or 
generated in a manner described below with reference to Figure 4. Alternatively, 
a candidate voxel set may include or be generated using all, essentially all, or 
nearly all voxels within (1) an entire imaging dataset; (2) a single or multiple image 
slices; and/or (3) one or more subsets or segmented portions of the imaging 
dataset. 

[0029] The hierarchical imaging information analysis procedure 300 may 

additionally include a VOI selection procedure 306 that involves selecting a first or 
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a next VOI for consideration. Additionally, depending upon embodiment details, 
the hierarchical analysis procedure 300 may include a washout analysis 
procedure 330, and possibly a plateau analysis procedure 320 and/or a persistent 
enhancement analysis procedure 310. 

A washout analysis procedure 330 may include a washout curve search or 
analysis procedure 332 that involves determining whether one or more candidate 
voxel sets within a VOI under consideration exhibit dynamic behavior that 
corresponds to a washout curve 130. In general, the washout curve search 
procedure 332 may determine whether one or more parameters corresponding to 
a candidate voxel set indicate an adequately, significantly, and/or strongly 
decreasing imaging signal intensity or relative enhancement profile over time. In 
certain embodiments, the washout curve search procedure 332 may involve (1) 
calculating or retrieving SER values for one or more candidate voxel sets within 
the VOI under consideration; and/or (2) determining which candidate voxel sets 
under consideration have SER values corresponding to or indicative of a washout 
curve 130. In one embodiment, a candidate voxel set may exhibit washout 
behavior or correspond to a washout curve 130 in the event that an SER value 
corresponding to the candidate voxel set is greater than a washout SER threshold, 
for example 1.1. A washout threshold may be specified, established, and/or 
determined in accordance with clinical conditions and/or statistical considerations. 

A candidate voxel set may additionally or alternatively be defined to exhibit 
washout behavior in the event that a slope and/or an angle corresponding to an 
imaging signal intensity curve or profile is negative or less than a washout slope 
and/or a washout angle threshold. For example, an exemplary washout slope 
threshold may correspond to a percentage decrease in imaging signal intensity of 
10% over a 5 minute interval, that is, a slope of -2 percentage units/min. 
Similarly, an exemplary washout angle threshold may be -5 degrees. | 

In the event that one or more candidate voxel sets currently under 
consideration correspond to a washout curve 130, the washout analysis procedure 
330 may initiate or execute an identification procedure 334 that involves 
identifying, selecting, and/or indicating as a malignancy or a likely malignancy 
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which of such candidate voxel sets exhibits a strongest, most significant, best, or 
most adequate correspondence to malignant behavior. Such correspondence to 
malignant behavior may be based upon a highest PE value, a strongest or most 
significant decrease in imaging signal intensity relative to a time period, and/or a 
most rapid decrease in imaging signal intensity relative to a time period. In one 
embodiment, an identification procedure 334 identifies as a likely malignancy a 
candidate voxel set corresponding to a highest PE value. Depending upon 
embodiment details, the identification procedure 334 may numerically, textually, 
chromatically, and/or graphically or visually identify a likely malignancy, for 
example, upon a display device. In particular embodiments, the identification 
procedure 334 may additionally or alternatively generate a graphical or visual 
representation of a washout curve 130 that corresponds to the voxel set identified 
as a likely malignancy. That is, the identification procedure 334 may generate, 
draw, or and/display a percent enhancement and/or signal intensity versus time 
curve based upon the imaging data that corresponds to the voxel set identified as 
a likely malignancy. Such a curve may be displayed separately or in an adjacent 
or overlapping manner relative to one or more portions of a corresponding 
displayed VOI. 

In one embodiment, an identification procedure 334 may identify more than 
one voxel set as a likely malignancy, possibly in accordance with an ordering or 
prioritization scheme. An ordering scheme may identify or indicate a voxel set 
having a highest likelihood of corresponding to a malignancy, a voxel set having a 
next highest likelihood of corresponding to a malignancy, and possibly one or 
more additional voxei sets in an analogous manner. In such an embodiment, an 
identification procedure 334 may visually specify or indicate each such voxel set in 
a unique manner, for example, using distinct colors. 

In the event that washout behavior is absent or apparently absent within the 
VOI currently under consideration, the hierarchical analysis procedure 300 may 
initiate or execute a plateau analysis procedure 320. In one embodiment, the 
plateau analysis procedure 320 includes a plateau curve search or analysis 
procedure 322 that involves determining whether one or more candidate voxel 



-11- 



sets within the VOI under consideration exhibit time dependent behavior 
corresponding to a plateau curve 120. In general, the plateau curve search 
procedure 322 may determine whether one or more parameters corresponding to 
a candidate voxel set indicate a constant, unchanging, or generally unchanging 
imaging signal intensity curve or profile over time. 

In certain embodiments, the plateau curve search procedure 322 may 
involve determining which candidate voxel sets under consideration have SER 
values corresponding to or indicative of a plateau curve 120. In one embodiment, 
a candidate voxel set may exhibit plateau behavior or correspond to a plateau 
curve 120 in the event that a corresponding SER value is less than a washout 
threshold, for example, 1.1, and greater than a persistent enhancement threshold, 
for example, 0.9 (i.e., the SER falls or generally falls within a range that indicates 
a flat or generally flat imaging signal intensity or relative enhancement profile over 
time). In a manner analogous to that for a washout threshold, a persistent 
enhancement threshold may be specified, established, or determined in 
accordance with clinical conditions and/or statistical considerations. 

A candidate voxel set may alternatively or additionally be defined to exhibit 
plateau behavior in the event that a slope and/or an angle corresponding to a time 
dependent imaging signal intensity curve or profile is less than a persistent 
enhancement slope and/or angle threshold (for example, +2 percentage 
units/minute and +5 degrees, respectively), and greater than a washout slope 
and/or angle threshold (for example, -2 percentage units/minute and -5 degrees, 
respectively). 

in the event that one or more candidate voxel sets currently under 
consideration correspond to a plateau curve 120, the plateau analysis procedure 
320 may initiate or execute an identification procedure 324 that involves 
identifying, selecting, and/or indicating as a malignancy or a likely malignancy 
which of such candidate voxel sets exhibits a strongest, most significant, best, or 
most adequate correspondence to malignant behavior. Such correspondence to 
malignant behavior may be based upon a highest PE value, and/or a most stable, 
unchanging, constant, or flattest imaging signal intensity relative to a time period. 
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In one embodiment, an identification procedure 334 identifies as a likely 
malignancy a candidate voxel set corresponding to a highest PE value. 
Depending upon embodiment details, the identification procedure 324 may 
indicate a likely malignancy in one or more manners described above. In one 
embodiment, an identification procedure 324 may identify more than one voxel set 
as a likely malignancy, possibly in accordance with an ordering or prioritization 
scheme, in a manner analogous to that described above. 

0038] In the event that washout and plateau behavior is absent or apparently 

absent within the VOI currently under consideration, the hierarchical analysis 
procedure 300 may initiate or execute a persistent enhancement search or 
analysis procedure 310. In one embodiment, the persistent enhancement search 
procedure 310 involves identifying, selecting, and/or indicating a candidate voxel 
set within or corresponding to the VOI currently under consideration that exhibits a 
strongest, most significant, best, or most adequate correspondence to malignant 
behavior. Such correspondence may be based upon a highest PE value, a most 
significant increase in imaging signal intensity relative to a time period, and/or a 
most rapid increase in imaging signal intensity relative to a time period. In certain 
embodiments, a persistent enhancement search procedure 310 identifies a 
candidate voxel set corresponding to a highest PE value as a malignancy or likely 
malignancy. The persistent enhancement search procedure 310 may indicate a 
likely malignancy in one or more manners previously described. In one 
embodiment, a persistent enhancement search procedure 310 may identify more 
than one voxel set as a likely malignancy, possibly in accordance with an ordering 
or prioritization scheme, in a manner analogous to that described above. 

[0039] The hierarchical analysis procedure 300 may further include a continuation 

procedure 350 that involves determination of whether another VOI requires 
consideration. If so, the hierarchical analysis procedure 300 may return to the 
VOI selection procedure 306; otherwise, further hierarchical analysis may not be 
required. 

[0040] Figure 4 is a flowchart of a VOI identification procedure 304 according to an 

embodiment of the invention. Depending upon embodiment details, the VOI 
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identification procedure 304 may include a screening condition determination 
procedure 402, a noise reduction procedure 404, a voxel set screening procedure 
406, and/or a volume connection procedure 408. 
[0041] A screening condition determination procedure 402 may involve 

determining, establishing, defining, identifying, generating, and/or retrieving one or 
more inclusion and/or exclusion thresholds that establish conditions that individual 
and/or combined voxels may or must satisfy to increase a likelihood of being 
diagnostically relevant. In one embodiment, a screening condition determination 
procedure 402 involves determining or defining a minimum enhancement 
threshold that individual and/or combined voxels must satisfy, for example, a 
minimum PE of 50%. 

[0042] Enhancement rate thresholds, such as achievement of a predetermined 

PE within a specified time period, may be used as a screening condition. A 
screening condition determination procedure 402 may additionally or alternatively 
involve calculating a difference threshold that may facilitate identification of 
individual and/or combined voxels having a valid or an invalid dynamic intensity 
profile. An exemplary invalid intensity profile may correspond to an individual 
and/or combined voxel for which S 0 = 2, Si = 4, and S 2 = 2. Although such a voxel 
exhibits 100% enhancement, it should be excluded from consideration due to its 
minimal or negligible maximum intensity (i.e., this is a very dim voxel). In one 
embodiment, a difference threshold may be determined by (1) calculating (Si - 
S 0 ) for each voxel to generate a subtraction data set; (2) computing a standard 
deviation for positive values within this subtraction data set; and (3) defining the 
difference threshold as a multiple (e.g., 3) of the subtraction data set's standard 
deviation. 

[0043] A noise reduction procedure 404 may involve generating combined voxels 

by averaging, blending, and/or combining imaging signal intensity information 
corresponding to a plurality of individual voxels to reduce or minimize the effect of 
signal noise. In one embodiment, the noise reduction procedure 404 involves 
averaging (or some other noise reduction technique such as Gaussian filtering) 
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precontrast, first postcontrast, and second postcontrast imaging slice information 
using a 3 x 3 x 1 or otherwise dimensioned kernel (e.g., 3x3x3). 

[0044] A voxel set screening procedure 406 may involve determining which 

individual voxels and/or combined voxels to retain and/or discard prior to a VOI 
selection procedure 306. In one embodiment, a voxel set screening procedure 
406 involves examining, comparing, and/or screening imaging signal intensity 
values corresponding to individual and/or combined voxels in accordance with 
inclusion and/or exclusion thresholds defined and/or determined in association 
with a screening condition determination procedure 402. A screening procedure 
406 may involve retaining those individual and/or combined voxels for which (1) an 
imaging signal intensity value exceeds a difference threshold; and/or (2) a PE 
value is greater than a minimum enhancement threshold. 

[0045] In one embodiment, a volume connection procedure 408 involves 

performing a connected component analysis upon individual and/or combined 
voxels that are retained as a result of a voxel set screening procedure 406, 
thereby generating one or more VOIs. A volume connection procedure 408 may 
additionally sort such VOIs in accordance with VOI size and/or other parameters. 

[0046] From the foregoing, it will be appreciated that specific embodiments of the 

invention have been described herein for purposes of illustration, but that various 
modifications may be made without deviating from the spirit and scope of the 
invention. For example, certain embodiments of the invention may determine 
whether a candidate voxel set exhibits persistent enhancement behavior in the 
event that washout behavior is absent, possibly by determining whether a signal 
enhancement ratio exceeds a persistent enhancement threshold and/or analyzing 
or evaluating imaging signal behavior relative to a persistent enhancement slope 
and/or angle. Accordingly, the invention is not limited except as by the appended 
claims. 
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